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We investigate by means of a number of different dynamical Monte Carlo simulation methods 
the self-assembly of equilibrium polymers in dilute, semidilute and concentrated solutions under 
good-solvent conditions. In our simulations, both linear chains and closed loops compete for the 
monomers, expanding on earlier work in which loop formation was disallowed. Our findings show 
that the conformational properties of the linear chains, as well as the shape of their size distribution 
function, are not altered by the formation of rings. Rings only seem to deplete material from the 
solution available to the linear chains. In agreement with scaling theory, the rings obey an algebraic 
size distribution, whereas the linear chains conform to a Schultz-Zimm type of distribution in dilute 
solution, and to an exponentional distribution in semidilute and concentrated solution. A diagram 
presenting different states of aggregation, including monomer-, ring- and chain-dominated regimes, 
is given. The relevance of our work in the context of experiment is discussed. 



Solutions of highlv elongated, cylindrical giant micelles are arguably among the best studied of the so-called equi- 
librium polymers [nl. Equilibrium polymers are formed in reversible polymerisation processes, and are in chemical 
equilibrium with each other, i.e., monomeric material is continually exchanged between the assemblies. An aspect 
not at all well understood is why ring closure seems to be unimportant in solutions of linear micelles, although this 
would remove unfavorable free ends ("end caps" ) from the solution. Closed loops have been observed in electron 
microscopic images of giant micelles |2j, but — as has been argued elsewhere (3|-|^ — in too low concentrations to 
significantly influence the properties of micellar systems. In other types of equilibrium polymer, such as liquid sulfur, 
the presence of rings is on the other hand thought to be all-important Q . It is believed that rings are suppressed in 
those self-assembled polymeric systems that are sufficiently rigid on the scale of the individual monomers . 

A useful model describing the self-assembly of giant micelles is what we call the Restricted Model of equilibrium 
polymers, where by definition ring closure and branching of chains are disallowed, and only linear chains form. In 
this model, the linear self-assembly is regulated by a free energy penalty associated with the free ends, the so-called 
end cap (free) energy E. The end cap energy is normally presumed to be a constant independent of the chain size or 
aggregation number, N , and of the overall monomer density (f) H. The basic scaling predictions for the equilibrium 



Despite their inherent polydispersity, self- assembled linear polymers resemble in many ways conventional polymers, 
that is, polymers with a fixed molecular weight. Indeed, the known statistical properties of conventional polymers have 
been applied quite successfully to predict the probability distribution function of the size of equilibrium polymers. 
In agreement with theoretical predictions 1^, the size distribution function has for instance been shown to decay 
essentially exponentially with size in the semidilute and concentrated regimes, because then correlations generated 



* emaibj wittmerQdpm . univ-lyon 1 . fr 



February 1, 2008 



PACS numbers: 82.35.-Ht, 61.25H, 64.60C 



I. INTRODUCTION. 




1 



by the excluded volume interaction are small |13|. A typical aggregate size distribution in a solution with strongly 
overlapping chains, obtained by computer simulation, is given in Fig. |l| (open squares). 

In the present work we relax the no-loop constraint, and discuss systems of (mainly) flexible equilibrium polymers 
under good-solvent conditions, where linear chains have to compete with rings for the available monomers at given 
total density 0, and end cap energy E. Branching of chains remains forbidden. We discuss in detail the number 
density of chains and rings of in different regimes (dilute versus strong overlap, ring versus chain dominated), and 
construct the complete "phase" diagram of this Unrestricted Model. Following up an an earlier, less extensive study 
[0, we report in this paper results obtained with a lattice-based Monte Carlo method |l^Jl^, and with a more recent 
off-lattice Monte Carlo scheme . We have succeeded in mapping the results of both methods onto each other, by 

making use of the natural length and energy scales that describe the configurations of the equilibrium polymers. This 
has enabled us to extract from the data the unknown prefactors that enter in the standard scaling theory, allowing us 
to construct a full diagram of states. Qualitatively, we confirm older, less elaborate descriptions y,^. As can be seen 
in Fig. |l|, the number density of rings as a function of aggregation number (indicated by the filled circles) is strongly 
singular, and dominated by a lower cut-off, that is, by the smallest ring allowed in the simulation. Effectively, this 
causes the crossover between dilute and semidilute regime, and that between ring- and chain-dominated regimes to 
coincide, at least for large end-cap energies E. 

This paper is organized as follows. We start in Sec. |l| with a brief presentation of the computational methods 
applied. Some technicalities concerning the parameters used and the configurations sampl ed ar e considered in Sec. III. 
The main computational results of this paper are presented in the Sections ^ and VII. In Sec. IV we discuss 
conformational properties, recall some notions and concepts of the physics of conventional pol ymer s (of fixed size), and 
determine the size ^ of the excluded volume blob . The following three Sections 0, ^ and VII are dedicated to the 
mass distributions. In the short Sec. ^ we fix some notions, refine the problem and pose specific tasks. Subsequently, 
we discuss the number density of linear chains of size N, ci(iV), in Sec. VI. Results obtained within the Restricted 
Model (with rings suppressed) and the Unrestricted Model (with rings allowed) are compared. We attempt to map 
both Monte Carlo methods onto each other, and discuss the physics at extremely high densities, where our off-lattice 
Monte Carlo approach exhibits features linked to packing effects, not present in the lattice-based Monte Carlo method. 
The number density of rings of aggregat ion n umber N, cq{N), and its connection to the size and mass distribution 
of the linear chains are discussed in Sec. VII - the core of the paper. In Sec. VIII we calculate the different regimes 
for systems of flexible equilibrium polymers in a good solvent, and compare this with the "measured" distribution 
of monomers over the rings and linear chains. In the final Sec. |^ we summarize our findings, and speculate on the 
relevance of ring formation in some recent experiments IfTHfT 



II. ALGORITHMS. 



All the computational algorithms we discuss are standard Monte Carlo schemes which have already been described 
and tested in depth in recent publications |p2[-p|. We shall therefore be extremely brief. 

The first method which was introduced to simulate the equilibrium polymers within the Restricted Model, is a 
grand canonical lattice algorithm based on a mapping on a Potts model ||T2|. Because this model lives on a simple 
cubic lattice, closed loops are by geometry forced to consist of an even number of monomers. Results obtained with 
this model will only be briefly mentioned (in connection with Fig. We discuss in more detail results obtained with 
two canonical Monte Carlo approaches — a lattice scheme [|l3| based on the Bond Fluctuation Model (BFM) [ p5[ , 
and a more recent off-lattice Monte Carlo (OLMC) approach |l4j generalizing an efficient bead-spring model. All 
BFM and OLMC simulations have been done in three spatial dimensions, the Potts model simulations in two and 
three dimensions. 

The BFM used in the present investigation is athermal, besides the constant scission energy Ubond = — J , which 
characterizes the bonded interaction: an energy J > is released every time a bond forms. Apart from the no- 
overlap conditions modelling hard-core excluded volume interactions, no other bonded (such as a stiffness potential) 
or non-bonded interaction has been used in the present study, though they might be readily included. 

In the off-lattice model, beads interact via a so-called (shifted) FENE potential for the bonded, and a purely 
repulsive Morse potential for the non-bonded interactions |T3] . The bonded interaction 



Ubondir) = -K{r„ 
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depends non-trivially on the distance r between two monomers. Here, J denotes again the constant part of the scission 
energy, and K the spring constant. The former quantity will be related with the model- independent effective end-cap 
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free energy E, in Sec.V]. Note further that the FENE potential is harmonic near its minimum at tq, while exactly 
at r = rg, Ubond = — J ■ The potential diverges logarithmically to infinity if r ^ Tmax and if r ^ Tmin = 2ro — rmax, 
where rmax and rmin are the maximum and minimum extension of the spring. Following ref. [ p^ , we set Tmax = 1, 
fmax — = To — Train = 0.3 and K — 40. Units are chosen such that fc^T = 1. 

Both models use local jump attempts. The time is measured, as usual, in Monte Carlo steps (MCS) per monomer. 
Every monomer is chosen at random, and allowed to perform a move subject to a Metropolis acceptance probability 
[^,0. In simulations the bonds between neighbors along the backbone of a chain are constantly subject to scission 
and recombination events. Since chains are only transient objects the data structure of the chains can only be based 
on the individual monomers, or, and this is the way we operate, on the notion of saturated and unsaturated bonds. 
An unsaturated bond does not connect a monomer to another one, a saturated bond does. Hence, we use no direct 
chain information, but lists of pointers linking up the bonds . With a given frequency one of the bonds (saturated 
or unsaturated) is chosen at random. If that bond happens to be saturated an attempt is made to break it, if it is 
unsaturated, i.e., if the monomer is at the end of a chain or a free monomer, an attempt is made to create a bond with 
another monomer sufficiently close. Obviously, for reasons of detailed balance, the bond formed must come out of the 
same set or range of bonds out of which also scission events are allowed to take place [|l3|,0 . Wc do not allow for 
branching in the present study. The minimal length of a closed loop is Nc = 3, i.e. closed dimers are not permitted. 
Free monomers are also not allowed to self-saturate their two unsaturated bonds. 



III. PARAMETERS AND CONFIGURATIONS. 



The only two model parameters of operational relevance in the present study to tune the system properties are the 
scission energy J and the number density </>. The starting configurations consist in both sets of simulations of randomly 
distributed and non-bonded monomers, which we cool down step by step (a sequence of so-called 'T- Jumps'), each 
step sampling a higher scission energy up to a maximum J = 15. This was done in order to produce a sufficient chain 
length and concentration variation to be able to put to the test the theoretical scaling predictions. 

Due to the constant breaking and recombining of the bonds, equilibration is much faster in our equilibrium polymeric 
system than is usually observed in systems of conventional polymers with fixed bonds. The algorithm presented above 
(using the pointer lists between bonds) allows us to simulate a large number of particles at very modest expenses of 
operational memory. 

In our BFM simulations, we varied the number density over three orders of magnitude from (f) ~ 0.000125 (1,000 
monomers per box) to </> = 0.075 (containing 75,000 monomers per box). For densities smaller than (j) — 0.0125 cubic 
lattices with volume 200"^ were used, whilst for higher densities a smaller box sufficed of volume lOO'^. Note that in 
the BFM every monomer consists of an elementary cube, whose eight sites on the cubic lattice are blocked for further 
occupation As a result of this, the volume fraction of material is given by 8^. As was shown elsewhere pO| , pl| , 
volume fractions of around 8</> « 0.5 are already quite dense within the BFM. The reason is that at higher densities 
than this, the system turns glassy due to the blockage of neighboring sites by other monomers. 

Most of the off-lattice Monte Carlo results involve 65, 536 particles for number densities between (f) = 0.0325 and 
(f> = 2 (and appropriately chosen box sizes) . Generally, we have sampled with the off-lattice Monte Carlo method more 
systems in the high and extremely high density regime, while we have explored with the BFM more systems in the 
dilute and semidilute regimes. Note that the highest densities in our off-lattice Monte Carlo simulations correspond to 
extraordinarily concentrated solutions. Indeed, if we estimate the corresponding volume fractions, we find that in our 
simulations these vary between 0.0072 and 0.44. (To obtain this estimate, an effective bead volume v « ttP /6 « 0.22 
was used, with a measured mean bond length of I « 0.75.) The latter value of 0.44 has to be compared with the (only 
slightly larger) hard-sphere freezing volume fraction ( "Alder transition" ) of about one half. Our simulations therefore 
extend to the "melt" regime of a dense liquid. 

In passing we note that strictly speaking the mean bond length / is not a constant, but decreases weakly with 
increasing density. However, this does not pose a serious problem, for we find that the total interaction energy per 
bond does remain roughly constant, with a Utotai ~ —(J— 0.4) for all overall monomer densities (f> and scission energies 
J probed. 

Periodically, the whole system is examined and various moments and distributions, such as the number densities 
Ci(iV), are counted and stored. Here, and below, we refer to ring-related quantities by using a subscript i — 0, and 
to chain-related ones by a subscript i = 1. The number densities Ci{N) are normalized such that = '^j^ Nci{N), 
with the overall density of monomers in species i. Obviously, the sum of monomers in rings, 0o, and that in linear 
chains, 0i, is equal to the overall monomer density (j) = 4>o + i'l- (Free monomers are counted as linear chains of 
length N = 1.) Obviously, (/) = (pi and 0o = in the Restricted Model, where rings are disallowed. We emphasize 
that the mean chain length (N) = n -^Ci(A^) / J^i n ^iiN) remains always two orders of magnitudes smaller than 
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the total particle number within the box. From our previous studies within the Restricted Model jl^jlj], we expect 
finite box-size effects to be small. 

Because of the differences in the lattice and off-lattice algorithms, we cannot expect the results obtained using 
these two methods to be directly comparable, even when the same parameters (J, 0) are used. To make a comparison 
possible, all relevant system parameters have to be mapped onto each other. This we do by taking the dilute limit as 
reference state. As the intrinsic energy scale we use E — J + 6J, with a model-dependent shift parameter 5 J presented 
in Sec. VI. The intrinsic length scales are the mean bond length I, and a length Ip to be discussed in the next section. 



IV. DISTRIBUTIONS OF SIZE. 



Let us first discuss the conformational properties of equilibrium polymers, and show that not only they follow the 
same universal laws as conventional polymers, but also that there is no essential difference in the behaviour of linear 
chains in the Restricted and Unrestricted Models. 

The configurational behaviour of the equilibrium polymers is most easily demonstrated by plotting their mean 
size as a function of the aggregation number N, as is done in Fig.^j. Indicated are the mean end-to-end distance 
of the linear chains, Rei{N), and the radii of gyration of the chains and the rings, Rgi{N) and Rgo{N). Averages 
have been taken over all linear chains {i = 1) or rings (i = 0) of given mass N in the simulation box. Only results 
obtained with the bond fluctuation method are shown, for two different densities and a fixed value of the scission 
energy J = 8. Similar results have been obtained with the off-lattice Monte Carlo method. Symbols are used for the 
distributions from the self-assembled chains in the presence of rings, lines for those without rings. Despite the fact 
that the simulations point at the presence of a large number of rings within the Unrestricted Model, with the dilute 
systems even being ring dominated, there is no measurable difference between the results for the linear chain sizes in 
the Restricted and Unrestricted Model calculations. This we quite generally find for all densities probed, and within 
both simulation methods. 

Indicated in the figure are two dashed lines, giving the theoretical slopes valid for very long polymers in dilute 
and concentrated solution In the dilute limit the chains are swollen and their size R is described by the scaling 
relation R (x N'^, with a self-avoiding walk exponent oi v = 0.588. Surprisingly, the relatively small rings shown in 
Fig.^ also closely follow the (asymptotic) scaling behaviour of their linear counterparts. In the following, the dilute 
limit will be used as reference state to be able to make a comparison between the different simulation models. As 
natural lengths we use the mean bond length I and an effective bond length b — bei, the latter determined from a fit 
of the simulation data to the scaling relation Rei{N) = beiN'^'^^^ ■ These two lengths in turn suggest a measure for 
the stiffness of the chains, which we shall call the "persistence length" Ip = b/l of the chains, although obviously it is 
a dimensionless quantity. For the BFM we find b ~ 3 and for the off-lattice Monte Carlo method b « 0.92. In 
both cases this gives a persistence length Ip = b/l that are of similar value (/ = 2.733, Ip ^ 1.2 for the BFM; / = 0.758, 
Ip « 1.1 for off-lattice model). The found values for b are essentially identical to what was obtained previously for 
monodisperse linear polymers ||l^,^. In a similar fashion, one may define and measure the prefactors bgQ and bgi 
from the radii of gyration RgQ and Rgi- In line with our previous work on conventional polymers studied in [pT[| , we 
find bgo/b « 0.3 and bgi/b w l/\/6. This result holds again for both simulation methods. 



Figure Q shows that, as expected |11|, the excluded volume correlations in our equilibrium systems are screened 
out for strongly overlapping chains. More precisely, the chains become Gaussian chains of blobs of size ^(0) = bgig'^, 
each blob containing g monomers, with 

g^gW^Gib'^y'^^'^''^ (2) 

and G an as yet unknown prefactor. The value of G may be fixed using the classical definition for the crossover density 
of a monodisperse solution: Att/3 £,^(j> — g. This definition gives for the prefactor G w 5.2, a value that turns out to be 
numerically consistent with alternative estimates obtained from fits to the scaling relations Rei{N) ~ b^ig'^ {N / gY^'^ 
or Rg,{N)^bgig-{N/gy/\ 

A further test of the configurational properties of the equilibrium polymers is provided by determining the depen- 
dence of the first moment of the size distributions on the mean aggregation number (Ni) = J2n ^Ci{N) / J2n Ci(A^)- 
Our results are presented in Fig.^, where we have plotted the radius of gyration divided by the blob size y = Rgi/£, = 
Rgi/(bgig'^), versus the mean aggregation number divided by the number of monomers per blob x = (Ni) /g{4>). In 
order to calculate these reduced quantities, we used eq.(||) together with our previous estimate for G = 5.2. The figure 
shows that the mean chain sizes of the living polymers again follow the same universal master curve as conventional 
polymers |l|]. Again we observe that the ring sizes, indicated by the full symbols, follow the same scaling law as the 
linear chains, which is actually rather surprising considering their relatively small size. The fact that the two slopes 
that indicate the scaling behaviour expected in the dilute and semidilute regimes cross nicely at {x,y) w (1,1), shows 
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that our estimate of G is actually rather accurate. Alternatively, we could have used this crossover point between the 
two regimes to define the prefactor G. 

In conclusion, the conformational properties of equilibrium chains are (within numerical accuracy) identical to those 
of polymers of fixed length. This is in contrast to earlier speculation in the literature, where it was surmised that 
rings might strongly influence the screening of excluded volume close to the crossover from the dilute to the semidilute 
regime We find that the universal functions are neither altered by the polydispersity, nor by the presence of rings. 
The blob size ^ is a function of total monomer density <j) only. 



V. MASS DISTRIBUTIONS OF EQUILIBRIUM POLYMERS. 



Next we turn our attention to the probabilities of finding aggregates of a certain aggregation number. For equilibrium 
polymers that are long compared to the effective bond length b, the main departure from conventional theory of flexible 
polymer solutions is that the reversibility of the self-assembly process ensures that the degrees of polymerisation 
are in thermal equilibrium. This means that the distributions co{N) and ci(N) are not fixed, but minimize the 
thermodynamic potential fl[co{N), ci(A'^)] of the system. It is natural to attempt a simplified theoretical description, 
using a Flory-Huggins type of mean-field approximation 

1 C30 

n[co{N),c^{N)] =Y.T. ^'(^) {^og{c,{N)l^) + fiN + f,) (3) 

1=0 N=l 

where we have written the thermodynamic potential as a sum over the different species i — for rings and i — 1 
for linear chains, and over all possible aggregation numbers N. The factor in the logarithm enters for dimensional 
reasons, where we set / equal to the mean bond length of the chains in the dilute regime; all energy units are measured 
in units of thermal energy, fcsT = 1. This ansatz for the thermodynamic potential is strongly motivated by its success 
describing the properties of equilibrium polymers within the Restricted Model, where ring formation is suppressed 
[^,0. The first term on the right is the usual translational entropy. The second term represents a Lagrange 
multiplier or chemical potential which fixes the total monomer density = 0o (a*) + 4'i (m) ■ contributions to the 
free energy which are extensive or linear in N are absorbed in this Lagrange multiplier. The as yet not specified terms 
fi describe the free energy contributions not extensive in the degree of polymerisation of the rings and linear chains. 
In general these may depend on the interactions between different chains and chain parts, and as a rule differ in the 
dilute, semidilute and melt regimes. We stress that in the prescription of eq. (H) terms (such as virial terms) which 
are not conjugate to co{N) or ci(A^) need not be made explicit, as they are absorbed in fi. Inspired by the results of 
the previous section, the central assumption of eq. (|^) is that the two distribution functions are only coupled via the 
chemical potential that makes sure that the total amount of monomers is conserved. 

In equilibrium, the distribution functions functionally minimise the thermodynamic potential, Sn/Sci — 0, giving 

l^co{N) = exp(-/o(iV, l^) _ f^N)H{N - N^p)) (4) 
l^Ci{N) = expi-Eilp) - fi{N, b^) - fiN) 

where we set /o + 1 = fo{N,b^<j),lp) and /i + 1 = E{lp) + fi{N,b^(j)) for convenience. The free energy associated 
with the linear chains /i is split into a model- and density-invariant end-cap free energy E — J -\- SJ (containing the 
operational scission energy J and a shift factor SJ discussed below), and a remaining part that somehow describes 
excluded volume correlations The Heaviside function H{x) enforces a smallest possible ring Nc{lp), in effect a 
lower cut-off. In actual systems this cut-off may depend on factors such as the detailed chemistry of the equilibrium 
polymers in hand, and on their bending stiffness not explicitly modelled here. 

Before we are able to comprehensively analyse our computer simulation results aided by eq. (|^), our task is (1) 
to complete the mapping of the simulation methods already started in the previous section, (2) to show that the 
Flory-Huggins ansatz is indeed justified, and that the ring and chain distributions decouple, (3) to identify the free 
energy terms /o and /i, (4) to show that the fi are only functions of N and the total density b^ip as indicated, and (5) 
to determine the equilibrium values of (po/ 4>, and {Ni) for a given system (J, (f). In the next two sections |VI| and VII 



we consider the distribution functions of the linear chains and rings, and analyse the free energy terms fi (tasks 3 and 



4). The mapping (task 1) will be completed in Sec. VIB, where we relate the natural (and in principle measurable) 



energy scale E with the operational parameter J. We show there that the distribution functions of linear chains in 
the Unrestricted Model can be computed from those in the Restricted Model, if the density of monomers in chains 



1 is given (task 2). In Sec. VIE we study the density crossover scaling fo r Ni ((j)i, 0, E). The relative distribution of 



monomers in rings and in linear chains (task 5) will be considered in Sec. VIII 
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VI. MASS DISTRIBUTIONS OF LINEAR CHAINS. 



In this section we focus attention on the size distributions of the Unear equihbrium polymers, as determined by 
computer simulation, within both the Restricted and Unrestricted Model settings. A typical example is given in 
Fig.|l|, where results are shown obtained with the off-lattice algorithm at high densities, that is, in the limit of strongly 
overlapping chains (0 — 1.5). The distribution functions for the Restricted and Unrestricted Models are both to high 
accuracy pure exponentials. Under conditions where linear chains dominate in the Unrestricted Model {(f>i ~ 0), we 
find that results from the Restricted and Unrestricted Models converge systematically. This is generally true in the 
limit of strongly overlapping equilibrium polymers, because then linear chains dominate the population of equilibrium 
polymers. (See also below.) At variance with recently published results of a simulation study |2^], we do not observe 
any sign of an algebraic singularity in the distribution function of the linear chains, even at the extremely high volume 
fraction of = 2. Results (not shown) from the BFM simulations confirm all trends observed with the off-lattice 
Monte Carlo treatment. 

A more critical test of our statement that rings only marginally influence the distribution functions of the linear 
chains is provided in Fig.^, we have plotted the distribution of linear chains ci{N) against the natural scaling variable 
X = N/ (Ni). Apparently, in the regime of strongly overlapping chains (main figure) the data collapse onto Ci{N) cx 
exp (— A^/ (A^i)) for a whole host of scission energies J. (Note that chain statistics deteriorate in the tail of the 
distribution.) In the dilute limit we find a similar data collapse, but with slightly different slope, pointing at a 
distribution of the type ci{N) cx exp (— 7A'^/ (Ni)), with 7 a constant identified below as a critical exponent. See 
the inset of figure ^, where results obtained from the BFM simulations covering the dilute limit are shown. Identical 
distributions have been found in both limits for systems without rings | [l3| , [l^ . 

How can our observations be understood? We have demonstrated in a previous publication that within the Re- 
stricted Model, the free energy term /i in the expression for the distribution function of linear chains in eq. (|^) can be 
calculated from the scaling theory of conventional polymers in a good solvent ||ll[]. Agreement with simulation data 
turned out to be remarkably good [|l3|,Q. We show here that t his an satz is also justified when rings are present, as 
long as the overall monomer density is not extremely high (Sec. VI D). 



A. The Mean-Field solvent. 



Starting with the (hypothetical) conditions where mean-field type of behaviour dictates the distribution functions, 
/i is a constant, /i may then be absorbed in E, i.e., in that case we may simply put fi = 0. From eq. (|^) we thus 
read off that within mean-field theory, the distribution function must be a pure exponential l'^ci{N) = e'xjp{— E ~ ^N) . 
It is easily shown that for this type of distribution, the equality n — 1/ {Ni) holds, albeit exactly only in the limit of 
large mean aggregation numbers (A^^i) S> 1. For the mean-chain length we obtain 

(A^i) = A?3Vi''exp(5i;) (5) 

with an amplitude A — exp(l/2), and mean-field exponents t] = 5 = 1/2. Note that (A^i) in the Unrestricted Model 
is a function of (jji — (j)i (0, E) , and only indirectly dependent on the control parameter 0. 



B. The dilute limit. 



Although largely unimportant in the semi-dilute and concentrated regimes, correlation effects do matter in the 
dilute regime, as we have in fact already seen in our discussion of the conformational properties in Sec. IV. Using the 
known statistical properties of self- avoiding walks, we infer that the free energy of a chain /i decreases logarithmically 
with the degree of polymerisation N 

h=-{l-l)\og{N). (6) 

where 7 « 1.158 is the susceptibility exponent of the n ^ vector model in 3 spatial dimensions Eq. (^) leads 
to a Schultz-Zimm distribution for the linear chains 

l^ci{N) = N'^^'^eyi^i-E ~ ^iN). (7) 

which indeed is born out by our simulation results, presented in the inset of Fig. ^ [^^ , since in the dilute limit we 
have ^ = 7/ (A^i) provided (A^i) > 1 |,|l|. 
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As already advertised in the previous section, we use the simulation results obtained in the dilute limit to map the 
different simulation models onto each other. The natural energy scale E = J + SJ is established by fitting the data 
for ci{N) to eq. (0), using the equality /i = 7/ (Ni). Having done this for the Restricted and Unrestricted Model 
simulation at various values of J and cj), we find SJ — 1.6 for the BFM and SJ = 1.7 for off-lattice model. We stress 
that SJ is model dependent, absorbing the different physical behaviour of the models on a microscopic scale (see 
Sec.|l|). It can only be coincidental that for the models wc used the SJ's turn out to be almost identical. Whatever 
the reason, the fixing of SJ completes the mapping. 

The concentration dependence of the mean aggregation number in the dilute regime, for future reference denoted 
by Nd, is given by 

Nd = {Ni) - j/fi = ^dZ^Vi" eMSE) (8) 

with the exponents = 5 = 1/(7+1), and a prefactor Ad — 7r'^(7 + 1) ~ 1.2, in terms of the exponent 7 and the 
usual gamma function F. Our simulations confirm these exponents and the prefactor (results not shown). 



C. Semidilute solutions. 

The aggregation number dependent free energy contribution eq. (|^) describes only dilute chains, i.e., chains which 
are too short to overlap. From the standard theory of conventional polymers jl^, one expects excluded volume effects 
to be screened out when the chains strongly overlap. Even then, this happens only for chains larger than the blob 
size, that is, for aggregation numbers N ^ g {<!>), where /i levels off to 

/i = -(7 - l)log(gi(0)) = /. + ^L_liog(630). (9) 

and gi{(j)) oc g{(f>)', the quantity gi scales like g (see eq. but with the different prefactor 

Gi =exp(-/,/(7-l))«0.43 (10) 

fixed by a constant fs — 0.13, estimated below. 

That /i indeed switches from length dependent (dilute) to density dependent (semidilute and melt) is shown 
in Fig. 1^, where we plot measured values of fi{N), using the mean- field relation l^ci{N) = exp{—E — fxN) and 
11=1/ (Ni). Figure | shows that /i is a constant of the degree of polymerisation at high overlap concentrations 
psf , as it should. It explains the scaling relation observed in Fig. ^, not least because the values of /i found for the 
Restricted and Unrestricted Models turn out to be identical if the overlap is strong enough. At low densities {(p < 0.5), 
/i becomes dependent on the chain length, due to the unscreened excluded volume interactions, but also on whether 
in the model ring formation is allowed or not. In the Unrestricted Model ring formation is allowed, and as a result of 
that some of the available material is stored in rings. This causes a shift in the distribution of linear chains relative 
to that in the Restricted Model. The trends of figure ^ confirm this. 

Using the measured plateau values of /i, we verify that the density variation predicted by eq. (0) holds, and 



estimate the prefactor Gi via eq. (10). This is done in Fig. |6|, where we have plotted the plateau values /i versus 
the dimensionless overall concentration of monomer, b^(j). Data from the BFM (pluses) and the OLMC (asterisks) are 
included in the figure. Also shown are /i values as they may be measured via the mean aggregation number in the 
strong overlap limit (SOL), that is, in semidilute and concentrated solution 

NsOL = {Ni) = ;3/20ii/2g^p f^^E+^h{b'q^)^ . (11) 

which depends explicitly on both and 4>i{<p,E). We use the directly measured (iVi) and (pi to obtain fi{b^(j)). 
This procedure gives identical results as the fi obtained straight from the distribution functions. Again, there is no 
observable difference between the results obtained with and without ring formation. In line with the prediction of 
eq. (^ ), we find a logarithmic density dependence (dashed line) with a prefactor (7 — 1)/(3j^ — 1) « 0.2. That in the 
semidilute regime both off-lattice and BFM results coincide reinforces our belief that the mapping between the models 
is robust. The divergence of t he lat tice and off-lattice data at very high densities is not problematic. A discussion of 



this issue we postpone to Sec. VI D , 

Within the semidilute regime, i.e., within the validity of eq. (^), the mean degree of polymerisation of the linear 
chains is power law function of cj) and (pi , and may be brought under the generic form of eq. ^ 

Ns = (Ni) = A,</.iZ3')0''-i eMSE) (12) 
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with (5 = 1/2, = (1 + (7 - - l))/2 « 0.6, and an amplitude 



3(7-l)/(3i/-l) 



1.1 



(13) 



with « 1.1. The (very weak) ^p-dependence arises because /i depends on 6'^^ rather than on l^(j>. Note that in 
the hmit of hnear chain dominance {<j>i/(f> ~ 1 ) one recovers the density dependence of the Restricted Model, i.e., 
Ni cx cjP'^ ||l,|l^- We verified the validity of eqs. ( |l^ and (|l3|) by comparison with the computer simulation data, but 
do not pause here to elaborate on the details. In Sec. VI E| we present the full crossover scaling of (A^i) , covering the 



entire range of concentrations from the dilute to the melt regime. 



D. Concentrated solutions. 

As we have discussed in a previous paper on the equilibrium polymerisation in the absence of ring formation [ p^ , 
eq. ([1^) only holds within the semidilute regime. This implies that there must be a large number of monomers g 
within a blob, for otherwise the blob concept becomes meaningless. Within our off-lattice Monte Carlo approach, we 
have probed such high densities that the semidilute description does break down - this is the melt regime already 
alluded to. It appears that we enter the melt regime if > 0.5, where a different physical behaviour intervenes, 
essentially due to fluid-like correlations resulting from local packing constraints. Clearly, effects of this nature cannot 
be expected to arise in any lattice-based bond fluctuation technique, because of the presence of an underlying lattice 
structure suppresses to a great extent fluid-like correlations. 

Our off-lattice Monte Carlo simulations point at a linear relationship between /i and the concentration material 
upon entering the melt regime: /i((/') = Bq + Bicl)b^ with Bq « —0.62 and Bi « 1.67. See Fig. ||, and in particular 
the inset. This also clearly shows that within the lattice model, the semidilute regime extends deeply into the melt 
regime, that is, even for concentrations where the blobs are so small as to contain no more than 5 « 10 monomers. 
We attribute this to the lack of a true fluid structure in the lattice-based model. 

We conclude from the above that while eq. (|^) breaks down in the melt, eq. ( pi] ) still holds (see Sec. VIE below). 



(A'^i) and (f> are then no longer related via a pure power law, but in addition via an exponential enhancement term: 
(iVi) cx 0^/^ exp . Direct measurement of (A^i) in our simulations conflrms this once more. Indications for 

the deviation from the usual power-law behaviour at very high densities have in fact also been found in earlier 
simulation studies ||2^,^ . The precise reason for the emergence of the exponential correction is at this point difficult 
to give. Theoretically, exponential corrections of the sort found here have been predicted theoretically for rigid and 
semiflexible micelles, but only within a second virial theory valid at low densities p7| . Within the second virial theory, 
the exponential density dependence of (A^i) arises from excluded volume interactions between the free ends and the 
central parts of the chains. Clearly, although the applicability of the concepts advanced in |^ to our high density 
system of flexible bead-spring chains is dubious, one may surmise that differences in the packing of the central beads 
and those near the ends of the chains could well lie at the root of problem. Further study is definitely warranted. 



E. Crossover scaling of the mean length of the linear chains. 

In this subsection we discuss the crossover scaling behaviour for the linear chain aggregation number (A^i), focusing 
on the (5-exponent (Fig. ^ and on the ?7-exponent (Fig. |^). We show that the computer simulation data from both 
Restricted and Unrestricted Models collapse onto the same universal function, if properly rescaled, and the directly 
measured density (pi is used. Again we include data from both the bond fiuctuation model (BFM) and the off-lattice 
Monte Carlo (OLMC) method. Results from the latter technique have been shifted upwards for reasons of clarity in 
both the figures |^ and ^. 

In Fig. 7 we compare the actual, measured (A^i) with a "hypothetical" dilute mean chain length /i(0i) calculated 
from eq. (8). In the dilute regime h = Nd, but outside this regime h represents an extrapolated value. Plotted is not 
(Ni) against h, but the quantity y ~ (Ni) /g^ against the quantity x = {g]^{(j))/h{(j)i)Y^^^^''\ specifically chosen to 
get a scaling with the 5-exponent: y oc x^ . Here, (/at is the mean chain length at the crossover density from the dilute 
to the strongly overlapping regime, defined by the equality g^ = N^, = Nsol- Equating eqs. (H) and (|ll|), gives for 
this quantity 

5^(0) = 4+2/(^-1) exp(-/i(<^)/(7 - 1)) (14) 

In the semidilute regime gN becomes proportional to g, with an amplitude Gn = A^^^^^''' exp [—fs/il — 1)] ~ 5 
very close to the blob prefactor G = 5.2. Because of the very large exponent 2/(7 — 1) « 12.7, we regard these two 
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values as numerically identical. Our conclusion is that gNi't') represents a generalization of the number of monomers 
per blob g{4>), also valid in the limit of very high densities where eq. does no longer apply. In passing we also note 
that Gn is also consistent with Gi (c/. eqs. (|l^) and (^). In other words, our characterisation of the excluded- volume 
blob is internally consistent. 

The purpose of all this is to make clear that in the high density regime discussed in the previous section, the mean 
aggregation number (iVi) can be described by the same universal function that is valid in the dilute and semidilute 
regimes, if the notion of the number of monomers per blob is generalised appropriately. The function qn that describes 
the generalised blob is universal in the sense that it depends solely on the overall monomer density (/>, and on the 
number of monomers per effective step length Ip. Indeed, the data points obtained for the Restricted and Unrestricted 
Models, within the bond fluctuation and off-lattice Monte Carlo methods, all collapse on the same universal function. 
The scaling plot confirms the scaling relation we set out to investigate y oc , with 5 — 0.46 in the dilute limit, and 
(5 = 1/2 in the semidilute and melt limits. 

In our second scaling plot. Fig. ^, we show the effects of a density variation on (A^i). To focus on the concentration 
behaviour of the critical exponent 77, we make use of eqs. (^) and (p^), and plot y — (Ni) jv^l'^ versus x — (pb'^v^/'^ 
where v = exp (E) 4'i/(j), and n — (3z/ — l)ip « 2.93 and <y9 = 1 -I- ( 1 -I- 7)/(3z^ — 1) « 3.81 are known critical exponents. 
This should yield a scaling relation y cx a;''. Equating equations (j8|) and (O) gives a crossover density I'^cf)* = Pv~^^''^ 
and a crossover length N* = Qv^l"^. The amplitudes P = [Aaj AsY^""^'^'^ « 1.9 and Q = A^P"^- w 1.6 are 
determined from Ad and As-, exactly as in the case of the Restricted Model, where v = exp(i?) |]l3|. The estimates 
cannot be expected to be very accurate because of the large exponent 1/(17^ — = 7.14, but are consistent with the 
estimates obtained by using the amplitudes Ad and Ag predicted from Gi (and/or G.) 

The data points for the dilute and semidilute regimes clearly collapse. As discussed earlier in the context of Fig||, 
the off-lattice Monte Carlo data do not conform to the scaling theory at very high densities, due to the effects of 
packing that dominate the melt regime. 



VII. MASS DISTRIBUTIONS OF CLOSED LOOPS. 



It has become clear from the discussion of Sec. [V, that the dimensions of the linear equilibrium polymers are 



successfully described in terms of the Flory exponent v. In the preceding Sec. VI we demonstrated that another 
critical exponent, the susceptibility exponent 7, is the exponent relevant to the description of the length distribution 
of the linear chains in the dilute and semidilute regimes. The length distribution of rings, only formed in the 
Unrestricted Model, is dominated by yet another critical exponent, a. This exponent is the well-known specific heat 
exponent, related to the Flory exponent v via the hyperscaling relation a = 2 — Dv in D spatial dimensions. As we 
shall make plausible below, the ring distribution is quite accurately described by the quasi singular scaling relation 
co{N) oc iV-(3-"). 

Let us return to main Fig. |l|. The distribution of rings indeed follows a power law distribution with an a = 0.5 in 
the dense limit of the off-lattice Monte Carlo simulation. This we expect, because in the dense limit = 1/2 and the 
simulation was done in I? = 3. The hyperscaling relation is checked extensively in the inset of the figure, where we 
plot the distributions as obtained for D ~ i not only by means of the off-lattice and bond fluctuation methods, but 
also from the grand-canonical Potts model simulations, mentioned in Sec. |ll[ Also included are data from the Potts 
model in two dimensions, where a = 1. Because the Flory exponent is = 1/2 for the concentrated or melt regime 
in both two and three dimensions, our results confirm the hyperscaling relation. 

The strong power law behavior seen in Fig. |l| does not exclude an additional exponential damping, which we expect 
to be present from eq. (^. In the strong overlap regime this exponential damping is difficult to detect, because there 
fi — 1/ (Ni) is very small. However, the exponential damping of the ring distribution is observable in the dilute 
regime, as is made clear by Fig. ^. In the figure we present data from the off-lattice simulations, for different s cissio n 



energies J and at a fixed concentration = 0.125. Note that (j^o increases with J, as does iVi = 7/^ (see Sec. VI] _ 
For large values of J we recover the power law behavior seen in Fig. l|, albeit with a different exponent a = 3 — 2.764 
because interactions are now not screened. This value is again in line with the hyperscaling relation, since in the 
dilute regime the Flory exponent obeys ly — 0.588. The drawn lines are fits eq. ([l6|). 

If the properties of rings and linear chains really decouple as was assumed in our Flory-Huggins ansatz, we can 
measure the free energy difference between rings and linear chains directly by plotting the ratio 



^ =exp(/o(7V,0) - (E + hiN,^))). (15) 

versus the aggregation number N. This is done in Fig. |l^, where we present off-lattice Monte Carlo data taken at a 
density of ^ = 1.5, and BFM data at (j> — 0.0625 (middle curve), the former shifted vertically for reasons of clarity. 
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Both densities are in the strong overlap regime. The lowest curve gives BFM data, taken at the dilute densities 
4> = 0.0125 (open symbols) and (j) = 0.00125 (filled symbols). The lines indicate power laws exp(/i — /o) oc iV^'^. 
From what we have said above about a, we expect t = 7 + Diy = 2.92 in the dilute limit, and t = 1 + D /2 = 2.5 
in the strong-overlap limit. The Monte Carlo data confirm this expectation, demonstrating once more the validity of 
the Flory-Huggins ansatz (task 2). 

Let us now briefly pause at a simple argument due to Porte, by which one may derive an expression for the ring 
distribution, and arrive at the exponent a [Q. The ratio of the ring and linear chain distribution functions, cq. ([T^), 
must be equal to the ratio of the respective partition functions, which in turn must proportional to the probability 
of opening a loop. The probability of opening a loop is proportional to (i) the Boltzmann weight exp(— i? — /i) to 
break a single bond |2^, (ii) the number of places where the ring can break, N, and (iii) the volume R^i that two 
neighboring segments can explore after being disconnected. Hence, 

^'co(iV)-Ao ^M~f^N) rj^j^^^y (Ig) 

V(i?,i(7v,0,g/-i)^ 



with Ao an unknown constant of proportionality. We have put eq. (|16|) to the test in Fig. 11, using the directly 



measured end-to-end distance Rei{N,(j)) discussed in Sec. IV, and the measured chemical potential /i, as discussed 



in Sec. VI. A wealth of data from both the lattice-based and off-lattice Monte Carlo methods in various regimes is 
included in the figure. The collapse of the data is next to perfect - the main result of this paper. The non-trivial 
behaviour of R^i {N, (j)) in the dilute and semidilute regimes explains the complex density dependence of the ring 
distribution in the various concentration regimes. The scaling plot yields a value of Aq ~ 0.1 for the constant of 
proportionality (Fig. ^), similar to the one estimated by Pfeuty and co-workers |6j in their analysis of self-assembled 
chains. 

While the scaling relation eq. ( p^ ) appears to be well satisfied for relatively large N, strong deviations from the 
universal asymptotic behavior is observed for small N. As seen from inspection of, e.g.. Fig. ^ or Fig. |o[ an 
unexpectedly large number of trimer rings are present in our simulations. This is emphasized in Fig. |l2|, where we 
give the fraction of monomers trimer rings 4>o{N = 3) relative to the total amount, and to that in all the rings, as 
function of the overall concentration of aggregating material. At densities below (j) — 0.25 most of all monomers are 
contained in trimers! The fraction of trimer rings seems largely independent of J. 

The effects of small chain length are systematically analyzed in Fig. ([l^), where we compare the measured free energy 
fa{N, (j)) with the asymptotic behavior /o"^™ — a~ [1 + Dv) log(A^), valid for large N . Here, a is a constant function 
of N, although it does depend on 0, because of the density dependence i?ei in the SOL. The density dependence 
shown in the inset, points at a decreasing tendency to form small rings with increasing concentration <p. It appears 
that the effect is much stronger in the off-lattice simulations than in the lattice-based approach, which reaches the 
asymptotic limit more rapidly (results not shown). 

It is important to emphasize that the observed small-ring behaviour is not in contradiction with the scaling picture 
presented in Fig. ^ It is in part caused by the similar deviation from the asymptotic behaviour of the dimensions 
of short linear chains, as can be seen in Fig. ^. We find that short chains are smaller in dimension than expected 
from the scaling behaviour for large N . As the entropy gain of opening a ring is smaller the shorter the chain, short 
rings become as a result of this more likely than expected. (See eq. (|6|).) In other words, the deviation from the 
asymptotic behavior of Cq {N) at small N is caused (at least in part) by the small- behaviour of the dimensions of 
the linear chains. 

That there are deviations from the universal asymptotic behavior in the small-ring population does not really come 
as a surprise. Unfortunately, even seemingly small shifts in the distribution of rings do matter when it comes to 
determining the relative amounts of rings and linear chains: the essentially algebraic distribution forces most of the 
ring mass to be concentrated in the smallest rings — exactly where the universal description based on scaling laws 
becomes inaccurate. This is the reason why obtaining a universal diagram of states based on theoretical argument, 
attempted in the next section, is fraught with difficulty, at least in principle. However, for the purpose of getting a 
qualitative picture of the aggregated states of equilibrium polymers, scaling theory provides a sufficiently accurate 
basis. 



VIII. DIAGRAM OF STATES. 

If we accept the theoretical distributions of eq. (^ at face value, and augment these with the values for the 
parameters 6, G, Gi, Ag and Aq obtained by fitting to the results of our computer simulations, we can calculate 
a diagram of states or "phase diagram" for our system of flexible equilibrium polymers p9| . For this purpose, we 
calculate, using eq. (||), the overall densities of monomers in rings and in linear chains, 0o and (jii, as well as the mean 
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chain lengths, (Nq) and (Ni), as a function of the total monomer density </>, and of the end cap energy E. For the 
mean end-to-end distance -Re, and the free energy correction /i to create an additional chain end, we simply put 

Re{N, h^4)) = minibN", bg^iN/gf'^) (17) 
/i(iV, = max(-(7 - 1) log(7V), -(7 - 1) log(gi)) 



with the monomer number per blob g{b^4>) = 5.2(6^0)"^/'^'^"-^' and gi{b'^4>) = O.O83(7(&'^0). Obviously, eq. (18) is 
only an approximation to the full universal functions. 

The set of equations to be solved require a numerical evaluation, essentially because the sums cannot all be evaluated 
analytically. For a given E we vary the Lagrange multiplier /u, and obtain from eq.(|^) the densities 4'i and, 
therefore, also A complication is that the distribution functions themselves depend explicitly on the density (j), 
at least in the semidilute regime. The total density is of course not known a priori, and has to be evaluated for 
any given value of ^. In principle, the only way out is a recursive iterative scheme. Fortunately, matters become 
simplified significantly if one starts the calculation in the dilute regime (at large /i) where both Re and /i are density 
independent. Recursive iteration may then be circumvented by slowly decreasing the value of /x, updating and 5(0), 
and then using a forward scheme to determine the next /i-value. 

We first consider the relative distribution of the monomers over the rings and the linear chains in flexible equilibrium 
polymers, where we set the lower ring cut-off at N^, = 3, conform our simulations. In Fig. |lj we compare the 
concentration dependence of the fraction of monomers in rings, 4>q/ <j), as obtained from our simulation studies, with 
the numerical results of the idealized model defined above. In the inset we furthermore give the concentration 
dependence of the mean chain lengths (iV), {Ni) and {Nq). Considering that no additional fit parameters were used, 
the agreement, although not perfect, is actually quite reasonable. The general trend is well described by the model 
calculation. 

The full diagram of states is presented in Fig. |l^. We distinguish three regimes in the {E, 0)-plane: (i) one where free 
monomers dominate at low densities and end-cap energies, (ii) one where rings dominate, at intermediate densities and 
high end-cap energies, and (iii) one where linear chain dominate, at high densities. The crossover to the ring-dominated 
regime is defined by the equality 0o = 0i = (/>/2. For energies larger than E = 6.7, we find two roots: the concentration 
(pri demarcating the crossover between ring and linear chain dominance, and the concentration separating the 
monomer- and ring-dominated regimes. The crossover concentration is defined as that concentration when half 
the linear chains is in monomeric form, and half in polymeric form of iV > 1. For comparison we have also indicated 
in Fig. ^ data found by means of computer simulation within the BFM. Configurations with ring dominance are 
denoted by circles, systems where linear chains dominate by squares. Apparently, the simulation data corroborate 
our theoretical phase diagram. 

In the figure |l5| we have in addition drawn the three dilute-semidilute crossover densities (/>* {E) of relevance to our 
discussion. The thin dashed line denotes the crossover density 1^(1)* (E) ~ 2 exp(— £'/3.8), valid within the Restricted 
Model (where rings are absent). The presence of rings alters qualitatively the crossover density, for in the Unrestricted 
Model it becomes independent of E at large end-cap energies. This bounds the semidilute regime to relatively high 
densities in the Unrestricted Model, as was in fact already remarked by Cates |^ . See Fig. ^ The semidilute regime 
cannot extend itself deeply into the ring regime, basically because the rings are too short. This is shown by the 
dashed-dotted line, which indicates the crossover density for the linear chains, as defined by equating g{4>l) with (Ni). 
(This definition of (pl corresponds to the one given in Sec. VI.) Note that 0.055, and that the chains start to 



overlap slightly below the 0ri-fine. The rings shorter than g remain swollen at densities well above (pri- This can be 
seen from the thick dashed line, the crossover density calculated by equating g{(j) = 0*) with the mean chain length 
(N) of all chains (chains and rings). This crossover density has an asymptotic value l^cj)* 0.29 in the limit of large 
end-cap energies. Around the crossover from the ring- to the linear chain-dominated regime, the mean chain length 
grows rapidly from {N) fa Nq Nc to {N) « iVi, because then most (but not all) of the additional monomers are 
now included in the long chains. 



IX. DISCUSSION AND CONCLUSIONS. 



Based on extensive computer simulations in two and three dimensions, using a number of different simulational 
techniques, we conclude that the configurational properties of the self- assembled linear chains are not altered if ring 
closure is allowed. The same can be said about the shape of the probability distribution of the linear chains. We 
find that the probability of finding linear chains of a certain length drops off exponentially with length, at least for 
concentrations where the chains strongly overlap. In dilute solution, where the chains do not overlap, the distribution 
function is of the Schultz-Zimm form. Previously, the same distribution functions were found theoretically, as well as 
by means of computer simulations, in equilibrium polymerising systems where ring closure was suppressed. Is ring 
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closure allowed, we find the resulting distribution of rings to be essentially algebraic, albeit that in dilute solution the 
algebraic distribution is dampened an additional exponential length dependence. 

It seems that the rings formed in equilibrium polymerisation merely deplete monomers from solution, making fewer 
of them available for absorption into linear chains. This is experimentally relevant, because it shifts the crossover to 
the semidilute regime to significantly higher densities. Interestingly, of the amount of material present in the rings, 
a very large fraction resides in the smallest rings allowed - in fact much more so than expected from the universal 
scaling relations which seem to be valid for large enough rings. In our case the smallest rings possible are determined 
by an arbitrary cut-off. In reality, the smallest possible ring is likely to be determined by the chemistry of the system 
in hand, making the cut-off a non-universal quantity. An important factor determining the minimum size of a ring 
may well be the rigidity of the bonds connecting the basic building blocks of the equilibrium polymers, and possibly 
also the configurational properties of the building blocks themselves. 

It has been argued in the literature that rings shorter than roughly a persistence length are highly unlikely to 
form [^j^ . If it takes many monomers to make one persistence length, ring closure is suppressed in favour of linear 
chains. This might be the reason why rings seem to be unimportant in giant micellar systems, because each individual 
surfactant molecule does not contribute a great deal to the length of an aggregate. For giant micelles the number of 
monomers required to grow an aggregate the size of a persistence length may be very large indeed, even when the 
aggregate is quite flexible. The same need not be true in other types of equilibrium polymerising system, and for 
these ring formation must be relevant. Indeed, it is well known that in liquid sulfur rings play an important role. We 
speculate that the short Sg rings, which seem to overwhelmingly dominate the ring regime of liquid sulfur, represent 
the aforementioned lower cut-off. 

The question arises why the problem of ring closure has not provoked a great deal of interest, for instance in relation 
to other types of equilibrium polymeric solutions. A plausible reason is that it appears to be difficult to experimentally 
distinguish between rings and linear chains, although new experimental techniques may change this in the future. 
Meijer and co-workers have recently been able to distinguish by means of NMR methods between material in tight 
rings and that in other chains in their system of hydrogen-bonded equilibrium polymers ||l^ . Unfortunately, these 
authors have not yet undertaken any systematic studies of the ring formation, for the moment barring a comparison 
with our computer simulation results. 

Direct measurements of entire size distributions are also extremely rare in the field of equilibrium polymeric systems. 
In fact, we know of only a single study that has been published in the literature. Greer and co-workers quite rece ntly 
presented size distributions of the living polymerisation of poly(a-methyl styrene) in the solvent tetrahydrofuran [[l9| . 
The distributions were obtained by means of size exclusion chromatography after termination of the polymerisation. 
The data very clearly show a distribution that is algebraic for small degrees of polymerisation, crossing over to an 
exponential distribution for large degrees of polymerisation, pointing at the presence of small rings. (The experiments 
cannot distinguish been rings and linear chains, so the distribution contains the sum of contributions from linear 
chains and rings, if present.) A mean-field power-law distribution proportional to iV~^'^ turns out to fit the data 
quite well at low degrees of polymerisation. 

We are somewhat puzzled by the good agreement with the (mean-field) prediction, for the equilibrium polymerisa- 
tion reaction is initiated by sodium naphthalene. This initiator forms bifunctionally reactive polymeric anions, which 
in principle do not allow for a ring closure reaction. However, there are several mechanisms by which rings may form 
regardless. One mechanism is that of Coulombic association of the charged end groups. This is possible because of the 
presumably strong binding or localisation of the cationic counterions to the anionic end groups, the solvent being very 
hostile to charged species. The resulting effectively dipolar chain ends may themselves form larger clusters (dimers, 
trimers etc.) for the same reason, allowing in principle for all kinds associations, including rings. Another mechanism 
could be the presence of small amounts of oxygen, which in the termination reaction lead to the formation of polymer 
radicals which could close up. Even if rings are indeed formed in the living polymerisation studied by Greer, we feel 
a direct comparison with our simulations is not warranted. The reason is that these rings must then be formed by 
processed competing with the on-going linear anionic polymerisation, requiring a different model altogether. 
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FIG. 1. Probability Distributions Functions of self-assembled chains in the Strong Overlap Limit (SOL). Main figure: 
Comparison of the computer simulation results for the Restricted Model (RM) and the Unrestricted Model (UM), obtained 
by the off-lattice Monte Carlo method (OLMC), described in the main text. The open squares represented the distribution of 
linear chains in the RM, the filled symbols the distribution of rings (circles) and that of linear chains (squares) within the UM. 
Data shown are from simulations at a monomer density of = 1.5, and scission energy of J = 7. At such a high density, the 
UM system is dominated by linear chains {<j}i/<l> ~ 0.88), and the distribtutions ci{N) of both models are virtually identical 
exponentials ci{N) oc exp{—N/Ni), where for convenience the notation A^i = (A''i) for the mean aggregation number of linear 
chains is used. For our choice of density and scission energy, {Ni^^^ ~ 120, Ni^^' « 111). The ring distribution is well 
described by the power law cq{N) oc ]^~^^~'^\ with a slope 3 — a = 2.5 (bold line). Inset: co(N) in a double logarithmic plot 
for the OLMC, BFM and Potts Model algorithm (vertical axes shifted for clarity). 
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FIG. 3. Scaling plot of the scaled radius of gyration Rg/£,(<j}) versus scaled aggregation number Ni/g{(j>), for rings (filled 
symbols) and linear chains (small symbols for RM, large open symbols for UM). Ni is here the mean aggregation number of 
the respective quantity presented. All data shown are BFM results. 
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FIG. 4. Normalized distributions for linear chains Ci{x) versus the scaled aggregation number x = N/Ni, where again 
A'^i denotes the mean aggregation number of the linear chains. Main figure: Data collapse on the exponential distribution 
ci{x) = exp(— a) for the OLMC data within the UM, a.t <j> — 1.5. Inset: Similar plot showing the BFM data in the dilute limit. 
As for the RM |13||l^ , we find ci{x) oc exp(— 7a;) with 7 « 1.16 a critial exponent. 
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FIG. 5. Density crossover of the quantity — log(ci (A'^)i"^) — E — N/Ni, obtained for J = 7 from the Restricted (open 
symbols) and Unrestricted Model (filled symbols) OLMC simulations. For conditions of strong chain overlap, this quantity 
becomes chain length independent and equal to the linear chain free energy /i(0) discussed in the main text. 
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FIG. 6. Scaled aggregation number of the linear chains A^'i versus the dimensionless concentration 06^ from OLMC simu- 
lations with the scission energy put at J = 7. Values from both the Restricted (open circles) and Unrestricted Model OLMC 
simulations (filled circles) are included. Also shown are the (asymptotic) values, obtained independently from Ci{N) for the 
largest J available for a given density for OLMC (asterisks) and BFM (pluses). At higher melt densities (f) > 0.5, we find for 
the OLMC data a nonalgebraic dependence of (j). This can be clearly seen in the inset where we fitted the data to a function 
/i(0) = Ba + Bitpb^, with Bo ~ -0.62 and Bi ^ 1.67. 
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FIG. 7. Scaling attempt for the average aggregation number of linear chains A^'i, to extract the S exponent discussed in the 
main text. The measured values for A^'i are compared with gN{4>) (defined in the text), and plotted versus {h/gN{4>))^^^^^'''' , 
where h is the (hypothetical) mean chain length for swollen equilibrium polymers without mier-actions between monomers of 
different chains. The linear chain monomer density (pi was measured explicitly. Restricted and Unrestricted Model simulations 
based on both OLMC and BFM methods are shown. The OLMC data are arbitrarily shifted upwards for reasons of clarity. The 
unshifted points all collapse on the same master curve! This confirms validity of the mapping between the different simulation 
models, and the universality of the scaling with the quantity gN{(j}), which we identified as the generalized number of monomers 
per blob containing the directly measured fi{4>) (see Fig.H). This relation is not altered by the presence of the rings. Hence, 
the free energy of a chain end is again shown to be a function of (j) = (jio + 4'i only. 
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FIG. 8. The scaled mean aggregation number of linear chains A^i versus the scaled density (j), chosen such as to get a scaling 
with the ?7 exponent. As explained in the text, A''i and (f> are rescaled with powers of = exTp{E)(f>i/(j>. In the Restricted Model 
(open symbols), and for large densities in the Unrestricted Model, v reduces to the density- independent affinity v = exp{E). 
Again, the OLMC data are shifted arbitrarily upwards for clarity. 
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FIG. 9. Ring distribution (symbols) obtained from the off-lattice Monte Carlo method for the Unrestricted Model in the 
dilute limit at a single concentration 4> = 0.125, but for different values of the scission energy J. Double-logarithmic plot. The 
fits to the theoretical curves discussed in the text (and indicated by drawn lines), confirm that the ring distribution co{N) 
depends on E only indirectly via n, i.e., via Ni. The plots confirm the hyperscaling relation a + Du = 2, and fix the prefactor 
of the scaling theory Aq w 0.1 in the dilute limit. 
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FIG. 11. Scaling plot of the quantity y = l^co{N) exp{nN) versus x = N{Rei{N)/l)^), using the directly measured 
end-to-end distrance Rei{N) of the linear chains for various systems for different algorithms and regimes, as indicated in the 
figure. The data collapse onto a master curve y = Xq/x with Ao ~ 0.1. 
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FIG. 12. Variation of the fraction material in trimers (f)o{N = 3)/0 with overall monomer density (p, as found in the OLMC 
simulations. Inset: variation of the fraction monomers in trimer rings, relative to that in all rings <j)o, versus the overall density 
of monomers. For all but the highest density, the mass is mainly concentrated in tiny rings. The J-variation is very weak. 
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FIG. 13. Small-ring effect for in the Unrestricted Model. Data were taken by means of OLMC simulations, in the strong 
overlap limit. Shown is the free energy difference fo{N, <j)) — fo"^"^™ versus the aggregation number N, where fo = — log(co) — /xiV 
and jg^^f™ represents the known asymptotic value for large N. Main figure: Results for the density cj) = 1.0, and various scission 
energies J. Inset: Similar to main figure, but now a fixed scission energy J = 7 and various densities 4>. 



26 



0.8 



0.6 



0.4 



0.2 





10 




• OLMC 
□ BFM J 
A BFM J 

- E-8+1 

- E=10+1 



10" 



10" 



10' 



FIG. 14. Comparison of the concentration dependence of the fraction of rings as obtained by Monte Carlo simulations, 
and that calculated from the theory, described in the main text. The curves (for Ao = 0.1, SJ = 1.6 and Ip = 1.1) contain 
no additional fit parameters! Main figure: Ratio 4'o/4' of monomers contained in rings. Data from BFM (open symbols) and 
OLMC simulations (filled symbols) are indicated. Inset: (A*'), TVo and A^'i for J = 10. Obviously, (N) « A^o for 4'o/4' ~ s^nd 
for large (j) we have A^^i oc with rj — 0.6. Note that the growth of (A'^) is intricate, and could in intermediate regimes lead 
to a power law with exponent ri > 0.6, in agreement with some (but not all) experimental findings U%- 
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FIG. 15. A diagram of states for the Unrestricted Model, with Ao = 0.1 and Ip — 1.1. Indicated are three regimes where 
monomers, rings and linear chains dominate the population of aggregates, as explained in the text. The ring regime does not 
extend to within the semidilute regime above the (N) = g{(f>)-lme. Also shown are for what values of E and l^if> we find in the 
BFM simulations configurations where cjxj/tf) > 1 (spheres), and those where 4>i/4> < 1 (squares). 
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